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The random field Ising model in three dimensions with Gaussian random fields is studied at zero 
' temperature for system sizes up to 60"^. For each realization of the normalized random fields, the 

strength of the random field, A and a uniform external, H is adjusted to find the finite-size critical 
■ point. The finite-size critical point is identified as the point in the J/-A plane where three degenerate 

ground states have the largest discontinuities in the magnetization. The discontinuities in the 
magnetization and bond energy between these ground states are used to calculate the magnetization 
and specific heat critical exponents and both exponents are found to be near zero. 



I. INTRODUCTION 

The random field Ising model (RFIM) is among the simplest statistical mechanical models with quenched disorder 
but is still not well understood. It is presumed to describe equilibrium phase transitions in physical systems such as 
fluids adsorbed in porous media and diluted antiferromagnets. However, comparisons between theoretical predictions 
and experiments have been inconclusive because of the difficulty of equilibrating the experimental systems. For 
the three-dimensional RFIM, it is known that there is an ordered phase for sufficiently low temperature and weak 
randomness ||; |j . The standard picture |^, ^, Q is that the phase transition is continuous and is controlled by a 
zero temperature fixed point with three scaling exponents. The zero temperature (strong disorder) fixed point implies 
that controlled renormalization group calculations cannot be carried out so that information about exponents has 
come from numerical simulations, series analysis ^ and real space 11 and other approximate renormalization 
" " "il, miTl and, in fact, it 



group calculations ||l^ . There have also been suggestions that the transition is first-order 
is difficult to determine whether the magnetization vanishes continuously or discontinuously at the transition because 
the value of the magnetic exponent, j3/v is very small. 

Monte Carlo simulations [1^ |l^, |l^, ^ of the RFIM suffer from long equilibration times and have been limited 
to small systems. The validity of obtaining critical exponents from small systems has been called into question by 
simulations showing that for 24^ systems even qualitative features such as the apparent order of the transition 
vary from realization to realization. The difficulties of long equilibration times and small system sizes for Monte Carlo 
simulations have prompted a number of studies of the zero temperature RFIM ||l^, |l|, . Ground states 

of the RFIM can be determined efficiently by mapping to the maximum flow problem and then using a polynomial 
time algorithm to solve the latter . The assumption of these studies is that the zero temperature transition is in 
the same universality class as the transition at nonzero temperature. In this paper we consider the zero temperature 
RFIM phase transition. 

Estimates of many of the critical exponents for the three-dimensional RFIM are converging but there is still a 
problem with the specific heat exponent, a. Monte Carlo simulations ||T^ and some zero temperature studies |2^, 
find a to be quite negative, for example Hartmann and Young ||2^ find a = —0.65. On the other hand, a recent zero 
temperature study by Middleton and Fisher pl| concludes that a is near zero. Some experimental measurements of 



the specific heat [g7| 29| show no divergence and can be interpreted as a near —1 while other measurements [^0|, |3l|] 
yield a near zero and the experimental picture remains controversial |2^, |3^, ^ . Large negative values for a are in 
disagreement with the modified hyperscaling relation, a = 2 — (d—9)v, that is a central feature of the zero temperature 
fixed point picture. The relatively well established results that v is in the range 1.1 — 1.4 and that 9 is very close to 
3/2 imply that a is not much less than zero. Very negative values of a are also inconsistent with the Rushbrooke 
relation a -I- 2/3 -I- 7 = 2 since 7 is believed to be close to 2. 

In this paper we study the zero temperature phase transition of the three-dimensional RFIM with Gaussian random 
fields. Our two primary goals are to provide evidence that the transition is, indeed, continuous and to measure the 
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specific heat exponent, a. The novel feature of our approach is that for each reaHzation of the normaHzed random 
fields we fine-tune both the strength of the random field and a uniform external field in order to bring the system 
to its finite-size "critical point" (we refer to it as a critical point in anticipation of our result that the transition 
is continuous). Machta, Newman and Chayes ijl^ implemented a similar idea in their Monte Carlo simulations. 
We identify the finite-size critical point as the point where three degenerate spin configurations have the largest 
jumps in the magnetization. Critical exponents are extracted from the finite-size scaling of the discontinuities in 
the magnetization and energy at the finite-size critical point. Our results support the view that the transition is 
continuous with a and /3 both near zero. 

II. THE RANDOM FIELD ISING MODEL AT ZERO TEMPERATURE 

The random field Ising model considered here is defined by the Hamiltonian: 

7Y = - J ^ SiSj - A^^hiSi - H^^Si (1) 

<i,j> i i 

where J is the coupling strength, A is the strength of the random field, H is the uniform external field, hi is the 
normalized random field at site i and Si is the Ising spin variable at site i. < i,j > indicates a sum over nearest 
neighbor pairs on a three dimensional cubic lattice of linear size L with periodic boundary conditions. We take J = 1 
and the random fields to be Gaussian distributed with zero mean and unit variance, 

P{h,)^^cM-f)- (2) 
The two quantities of primary quantities that we measure are the magetization, m, 

m = -1 ^ (3) 

i 

and the bond energy, e, 

e = I] s^SJ. (4) 

<hj> 

The presumed phase diagram of the three-dimensional RFIM is shown in Fig. |l|. The solid line is the phase transition 
between the ordered and disordered phases. The point (Tc, 0) is the critical point of the pure Ising model while the 
zero temperature phase transition is at the point (0, Ac). Assuming the absence of special points along the phase 
transition line (e.g. a tricritical point) the universal properties along the entire phase transition line, except at A = 0, 
are expected to be the same as at the zero temperature transition. 

Consider the zero temperature transition. If A < Ac the system is in one of two ordered phases so the magnetization 
as a function of H has a jump at if = while for A > Ac, the magnetization is a continuous function of H. If the 
transition is continuous, the spontaneous magnetization, ip is expected to vanish as a power law as A approaches Ac 
from below, 

1^ ~ (Ac - A)0 (5) 

where ip = m{H 0+) and j3 is the magnetization exponent. Figure |^a illustrates the continuous transition scenario 
in the zero temperature H-A plane with a critical point at the end of a line of first-order transitions. Another 
possibility is that the zero temperature transition is first-order. A possible scenario is illustrated with Fig. Here 
[H = 0, A = Ac) is a point of coexistence of two ordered phases and one disordered phase. The magnetization of the 
coexisting ordered phases is nonzero while the disordered phase has zero magnetization. There will also be a nonzero 
"latent heat" at the transition. Although entropy is ill-defined at zero temperature, it is reasonable to define latent 
heat in terms of a discontinuity in the bond energy, e, between the ordered and disordered phases. 

The foregoing applies to infinite systems. The ground states of a typical finite system for a given realization of the 
random field are shown in Fig. ||. Each point in the H-A plane corresponds to a single ground state of the system. 
The set of points corresponding to a single ground state form a polygon since, for a given spin configuration, the 
energy is linear in both H and A. In principle, there might be several degenerate ground states in a nonzero area of 
the H-A plane but, for a continuous distribution of random fields, the probability of exact degeneracy vanishes except 
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along lines and at points. Along the edges between polygons the two ground states corresponding to each polygon are 
degenerate. Three ground states are degenerate at "triple points" where three edges meet. In principle, more than 
three ground states can be degenerate at a point but, for a continuous distribution of random fields, the probability 
of more than three degenerate ground states vanishes. 

Degenerate ground states typically differ on a small fraction of spins but, corresponding to the first-order line of the 
infinite system, some degenerate ground states differ by a large fraction of the total number of spins. The bold lines 
in Fig. ^ correspond to large jumps in the magnetization. For A < Ac, there is a single jump between ground states 
with large positive and negative magnetization and the separation between these "phases" is the piecewise linear 
curve 7?coex(A). This coexistence line is close to but not coincident with the A axis. The spontaneous magnetization, 
■0 is the magnetization of the positively magnetized ground states along i?cocx- As A^, is approached from below, tp 
decreases in steps at triple points. Since the net change in magnetization around a triple point is zero, the decrease in 
magnetization along -ffcoox is exactly the jump across the third edge or "lateral line" that defines the triple point. As 
the critical point is approached, the magnetization drops more precipitously so that the magnetization jump across 
the lateral line increases. We identify the finite-size "critical point" as the triple point with the largest jump across 
the lateral line. In Fig. |^ the finite-size critical point is shown as an open circle and the corresponding lateral line is 
also shown as thick line. The three states that are degenerate at the critical point are called the — and states 
according to their relative magnetizations. If the transition were first-order these states would be the zero temperature 
analog of -I- and — ordered phases and the disordered phase that coexist at a thermal first-order transition. 

One way to measure the discontinuity in the magnetization at the critical point is via the quantity, m*, 

m* — 2^^^^ ~ ™o)(™- ~ '71o)(to+ — m_) (6) 

where m+, toq and m_ are the magnetizations of the three coexisting critical ground states. Note that the same 
combination of three magnetizations is small at triple points away from the critical point, even along the first-order 
line. 

Although the energy, H, is a continuous function of H and A, the bond energy, e, has discontinuities across edges 
separating coexisting ground states. The bond energy along the coexistence line, i/cocx(A), is expected to increase 
monotonically as A increases with jumps across each lateral line. The largest discontinuity in the bond energy is 
expected at the finite-size critical point and the quantity e*, is a measure of that discontinuity, 

e* = i(e+ -t- e_)(eo - e+)(eo - e_) (7) 

where e+, e_ and cq are the bond energies of the — and states, respectively. 

If the phase transition is continuous, both m* and e* must approach zero as the system size goes to infinity while 
if the transition is first-order, these quantities will saturate at nonzero values. Furthermore, if the transition is 
continuous, we propose the following finite-size scaling behavior for the disorder averages of these quantities, 

^ ^ (8) 

and 

^ L-^Ci-")/'' (9) 

where L is the linear size of the system, a is the specific heat exponent, /3 is the magnetization exponent and v is the 
correlation length critical exponent. The finite-size scaling hypothesis for the magnetization is essentially identical 
to the standard finite-size scaling hypothesis except that the measurement is made at a point that is fine-tuned for 
the given realization of disorder rather than at the infinite system size critical point. The finite-size scaling for the 
bond energy energy is more difficult to motivate. First, at zero temperature free energy and energy are identical and 
entropy and specific heat are not well-defined. It is reasonable to suppose that the derivatives of energy with respect 
to J at zero temperature have the same singularities as the derivatives of the free energy with respect to temperature 
at nonzero temperature. Thus, from Eq. (Q), bond energy plays the role of entropy and its derivative with respect to 
J plays the role of specific heat. The conventional finite-size scaling hypothesis for the specific heat is 

C[T) ^ L-^/^CiiT -T^)L^/'') (10) 

If the specific heat is integrated across the finite-size rounding region defined by \T - < L^^/'" we obtain the 
finite-size "latent heat" , I of the transition. 
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At zero temperature, this latent heat is replaced by the discontinuity in the the bond energy. Furthermore, since 
the bond energy changes at a discrete set of discontinuities, we propose that a finite fraction of the latent heat is 
concentrated at the triple point with the largest discontinuity, i.e. the point we have identified as the finite-size critical 
point. Since e* measures the square of the latent heat at the finite-size critical point, we obtain Eq. fl). 



III. SIMULATION METHOD 



The main goal of our simulations is to find and measure the properties of the zero temperature, finite-size critical 
point. The finite-size critical point is defined as the triple point with the largest jump in the magnetization between 
all three coexisting states. To find it we carry out an iterative search using a method similar to the one described 
in Ref. |Q. First, a point in the ordered phase on the first-order line is located. For example, in Fig. || we locate 
ground states "a" and "b" and a point on the coexistence line between them. Next, the first-order line is followed 
in the direction of increasing A from one triple point to the next. In the example of Fig. ^ the triple point where 
"a", "b "and are degenerate is located first. At each triple point, the continuation of the first-order line (in this 
example, between "a" and "-)-" ) and the lateral line (between "b" and ) are identified according to the relative 
sizes of the discontinuities in the magnetization across the lines. The sequence of triple points is recorded and the 
triple point with the largest discontinuity in the magnetization across the lateral line is identified as the finite-size 
critical point. 

The most difficult computational problem in carrying out this program is to find the ground state for given values 
of A, _ff and {hi}. Finding ground states can be reduced to the problem of finding the maximum fiow on a graph 
p2| , |35| for which there exists polynomial time algorithms. We use the push-relabel network flow algorithm |36|, p7| , 
implemented in as version hi_pr, to find ground states. 

Let's now describe in detail the two subroutines of the search algorithm. Subroutine 1 finds a point on the first-order 
line separating the positively and negatively ordered phases. Two ground states are found deep in the ordered phases 
at values Aq = 2.2 and H = ±0.1. Once these ground states are found, the algorithm calculates the value of H where 
these two states are degenerate along the A = Ap line. This point of degeneracy is easily found by the solution of 
two simultaneous linear equations since once Ti., e and m are known for a given spin configuration at some point in 
the H — A plane, the value of Ti. for that spin configuration at any other point is linearly related. At this point of 
degeneracy the ground state is calculated. If the result is either of the two original ground states then we have found 
a point on the first-order line. If the ground state at the point of degeneracy is not one of the original states then 
a new point of degeneracy is found between this new state and one of the original two states. The choice of which 
original state to pick is made by the criterion that the jump in the magnetization between it and the new state is 
biggest. This process is repeated until the first-order line is found. 

Subroutine 2 follows the first-order line found in subroutine 1 until the first triple point is found in the direction 
of increasing A. From the point found by subroutine 1 on the first-order line a new point is picked by moving in the 
direction of the first-order line and increasing A by 0.1. On Fig. |^ this new point is on an imagined continuation of 
the thick line between "a" and "b" in the direction of An increase of 0.1 in A has proved to be sufficiently large to 
extend to a new ground state for all system sizes studied. Subroutine 2 is similar to subroutine 1 except that original 
ground states are always retained in the iteration and a sequence of new ground states with decreasing values of A 
are generated. At each step in the iteration, the point of degeneracy is found between the new and the old states and 
the ground state is computed at that point. If the ground state is one of already identified states then the triple point 
has been found otherwise the iteration is repeat. Call this first triple point {Hi, Ai). 

The triple point (iJi,Ai) defines three lines of coexistence, one of these lines is the previously identified first- 
order line. The jumps in magnetization across the other two lines are measured. The line with the smaller jump in 
magnetization is called the "lateral line" and the remaining line with the large jump in magnetization is continuation 
of the first-order line. The jump across the lateral line is recorded and subroutine 2 is invoked again to find the next 
triple point along the first-order line. By repeating this set of steps many times, a sequence of triple points {Hi, Ai) 
are identified. The triple point with the largest jump across the lateral line is identified as the finite-size critical point, 
{H*,A*). 

We note that while the critical point is identified by the size of the jumps in the magnetization, in every case that 
we examined it also has the largest discontinuity in the bond energy. In principle, a definition of the finite-size critical 
point based on the bond energy discontinuity might sometimes yield a different triple point. 

A typical running time for finding the ground state for a system of L = 10 is 6 x 10~^ sec/spin on a Pentium 
III 750MIIz machine. The total running time of the entire algorithm for the same system size is approximately one 
minute. We simulated systems of size up to L = 60. The number of realizations of disorder range from 23 to 126 for 
different system sizes. The small number of realizations are a consequence of the fact that many ground states must 
be explored for each realization of disorder to find the finite-size critical point. 
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IV. RESULTS 



Figure Hshows the critical magnetization discontinuity m* as a function of system size L. Clearly this plot, together 
with Eq. (H) suggests that the magnetization discontinuity at the critical point does not decrease with L which can 
be interpreted either as a first-order transition or a value of (3 very near zero. The quality of the data is not sufficient 
to make a useful measurement of (i. Figure ^ shows the finite-size scaling of the critical bond energy discontinuity e*. 
A fit of the form: 



e'^a + bL-" (12) 

yields a — 0.007 ± 0.005, 6 = 21 ± 5 and c = 1.7 ± 0.1. It is clear that e* vanishes as power of L. This result is a clear 
indication of a continuous transition. Assuming a=0, we obtain a fit, shown in Fig. ^ of the form: 

e* = bL-" (13) 

with b = 16.0 ± 2.6 and c = 1.59 ± 0.05. From the latter result and Eq. (||) we obtain (1 - a)/v = 0.80 ± 0.03 
(with — 2.84, jd.o.] . = 0.57 and Q = 0.73). This result is in agreement with Middleton and Fisher's value 
0.82 ±0.02. 

The correlation length exponent v and the infinite size critical disorder strength, Ac can be obtained from the finite 
size scaling of A*. The fit in Fig. ^ is: 

A* = Ac + bL-^'" (14) 

with Ac = 2.29 ± 0.02, & = 4.1 ± 0.6 and = 1.1 ± 0.1 (with = 2.79, x^/d.o.f. = 0.69 and Q = 0.59). This result 
is in relatively good agreement with recent results in the literature as shown on Table |. The value of v is somewhat 
smaller than the recent values in Refs. |2^, |l] which may result from the use of a smaller maximum system size in 
our study. Combining our results for v and for (1 — a)/i^ yields a = 0.1 ± 0.1. Using the larger values of v from 



Refs. 20 



21 gives slightly negative values of a. In any case, within the uncertainties, our results are consistent with 



modified hyperscaling. 



V. DISCUSSION 



The main result of this study is a new measurement of the specific heat exponent of the RFIM. Our result is a 
value of a near zero that is consistent with the modified hyperscaling relation and the recent simulation results by 
Middleton and Fisher pl| . It is not clear why a very similar measurement by Hartmann and Young pQ] give a quite 



negative. Our value of a and that of Ref. |£] 
critical point whereas in Ref. pQ] the average 



are based directly on the finite-size scaling of the bond energy at the 
3ond energy is numerically differentiated to obtain the specific heat and 



the finite-size scaling of the peak is used to obtain a. 
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Ref. 

This work 

|io| 

|L9| 

\m 



2.29(2) 
2.28(1) 
2.270(4) 
2.26(1) 
2.29(4) 
2.37(5) 



1.1(1) 
1.36(1) 
1.37(9) 
1.22(6) 
1.19(8) 

1-0(1) 



0.80(3) 
1.20 
0.82(2) 



0.12 
-0.63(7) 
-0.12 



1.55 -0.55(20) 



TABLE I: A summary of recent zero temperature estimates of Ac, i^, (1 — a)/]/, and a. For a and (1 — a)/i^ the value without 
error estimated is derived from the other, directly measured, value and the same authors' value of v. 
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FIG. 1; Phase diagram of the random field Ising model. The ordered ferromagnetic phase is labeled F and the disordered 
paramagnetic phase is labeled P. The curve is a line of phase transitions. 
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FIG. 2: Possible phase diagrams for the zero temperature phase transition: (a) continuous transition, (b) first-order transition. 
+ and — are the coexisting ordered phases and is the coexisting disordered phase in (b). 



FIG. 3: Ground states for a given realization of normalized random fields, {hi}. Coexistence lines between ground states with 
very different magnetizations are shown as thick lines on the plot, the finite-size critical point is shown as an open circle and 
the three coexisting states at the finite-size critical point are labeled -|-, — and as described in the text. 
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FIG. 4: The magnetization discontinuity at the critical point, m* vs. system size, L. 
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FIG. 5: Log-log plot of bond energy discontinuity at the critical point, e* vs. the system size, L. The solid line is a fit as 
described in the text. 




FIG. 6: The critical strength of randomness. A* vs. system size, L. The solid line is a fit as described in the text. 



